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SUMMARY 


Heat transfer data measured along th ® an^ STS- ^ ar ^ comp ^ ed"" ti t h 
fuselage of the Space Shuttle orbiter during STS-2 and bib J are <- t* 

predictions of empirical heating techniques derived from co f U ght 

Steps required to extrapolate an existing leeward centerline theory to flight 

rfppirimatelf oreseSfed! The side fuselage^ 

ST e r«^^ 

of embedded eoetlce. e^neg ^ tbe ■ « 
impingement line which significantly enhance local heating ra^es a 
tunnel and flight conditions. 


INTRODUCTION 


-Heating on top of the Space Shuttle orbiter 's vortex-dominated fuselage 

1, , complex function of Hech number . In t.ri'.f 
The upper fuselage thermal environment is generally character 

orbiter haf beeS developed an'd successfully applied to wind-tunnel^data ^ 
covering a large range in Refolds n^r ^ngle ^ attac. correlation 

10. 4 This method consists of a modified turbulent swept cy oil-flow 

using an effective loci .veep angle “ « Snntcted 

patterns on the upper fuselage. cons anzles and the distribution 

between the axial distribution of measured sweep angles and tne a 

of top centerline heating. Thin report ^im, number. 

wind-tunnel sweep angles to account for condl heating predictions with 

and Mach numbers. Comparisons of leeward centerline heating prea 

flight values are then presented. 

The basic concepts for a new technique designed to predi l5.g e ^^hod uses 
the side fuselage impingement line "* .‘’centeriire'theory. Furthermore, 
U STof SScS ‘t“.r conceding the relationship ^™en .«r- 

direction, and the aide fuselage Impingement he.tl-8 distribution. 
The side fuselage method is derived from oil-flow and phase-cha g P 
wind-tunnel data and supplemented by thermocouple measurements. Although 
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si/Je impingement method is still under development, some promising preliminary 
comparisons have been made with both wind-tunnel and flight heating rates. 


SYMBOLS 


AFFDL Air Force Flight Development Laboratory 

G Chapman-Rubes in coefficient 

L characteristic length of wind-tunnel model or full-scale vehicle as 

indicated 

M Mach number 

9 - 

CMS Orbital Maneuvering System 

<Ic convective heating rate 

^c h stagnation heating rate on a sphere with radius equal to average 

height of orbiter flat side body 

3c 0 stagnation heating rate on a scaled one-foot radium sphere in 

the free stream 

^Cg stagnation heating rate on a sphere with radius equal to that of 

orbiter top fuselage 

Re Reynolds number, based on L unless otherwise specified 

s cross-sectional surface running length measured from top centerline 

STS Space Transportation System 

T temperature ^ r 

x axial length measured from orbiter' s nose 

L centerline 

a angle of attack 

ot r reference a defined in equation (6) 

0 bov shock angle measured with respect to free stream direction 

5 flow deflection angle across bow shock 

: surface flow angle and local angle of attack 

1 change in sweep angle with respect to angle of attack 

i T change in n with respect to rj L 
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A effective sweep angl'' 

A r reference A defined in quations (3), (4), and (5) 

v Prandtl-Meyer turning angle 

$ Mach number correction factor 

X hypersonic viscous -interaction parameter = H 3 y' C/ y Re (see 

reference 5) 

Subscripts 


D 


ext 

i 

L 

l 


max 

t 


in equation (1), tv*ice the orbiter upper fuselage radius; in 
equation (16) , twice the average height of orbiter flat side body 

extrapolated 

flight condition 

quantity based on characteristic length 
quantity based on local flow properties 
maximum value 
wind-tunnel condition 
free stream 


WIND-TUNNEL DATA 


Oil flow patterns used to extrapolate the leeward centerline heating 
method to flight conditions were obtained on the upper fuselage of a 
0.01-scale orbiter model in air in the Langley Research Center’s Mach 6 and 
Mach 10 facilities , 6 ~ 8 and also on a 0.006-scale model in 20- Inch Mach 
14 AFFDL wind tunnel at Wright-Patterson Air Force Base. 9 Tests at Mach 6 
ranged from 15° to 40° angle of attack at Reynolds numbers of 2.7 x 10 & , 5.4 x 
10 & , and 7.3 x 10 6 . Oil flows at Mach 10 were run at a * 15° to 45° for free- 
stream Reynolds numbers of 0.59 x 10 6 , 1.19 x 10 & , and 2.37 x 10 & . Angles of 
attack of 15° to 40° at Re^ =* 0.280 x 10^ and 0.423 x 10 6 w< re used at Mach 
14. 


The phase-change paint heat transfer, oil flow, and thermocouple measure- 
ments used to derive the side fuselage impingement heating method were made on 
0.01-scale Shuttle orbiter models in Langley’s Mach 6 and Mach 10 tunnels. 
Tests were conducted for Re« ■ 0. 59 - 7.3 x 10 6 and at angles of attack of 
20°, 30°, and 40°. Oil-flow tests were made using an aluminum model, ftodels 
for the phase-change paint heating measurements were constructed of a filled 
epoxy casting compound and a semi-infinite °lab solution 10 was assumed during 
act a reduction. The supplemental thermocouple results were drawn from a 
previously unpublished data base described irl reference 4. 
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FLIGHT DATA 


Flight measurements used for comparison with the top centerline and side 
fuselage heating methods were obtained on STS-2 and STS— 3 at the locations 
shown in figure 1. The top centerline heating rate* . measured by calorime- 
ters, were the only data from STS-2 that were used In this report. The con- 
vective component of heat transfer for the calorimeters was determined by add- 
ing the radiative loss term to the calorimeter value. Keating due to solar 
radiation was then subtracted for those trajectory points where the 
instruments were In sunlight. The effect of solar heating was computed by the 
technique of reference 11. All of the instrument locations for STS-3 were ^ 
occupied by thermocouples. A one-dimensional, transient. -conduction analysis 
was used to determine convective heating for these Instruments with solar 
radiation heating, once again, computed separately.* A process combining the 
results of trajectory and atmospheric reconstruction provided information on 
vehicle attitude and free stream conditions. ’ 


LEEWARD CENTERLINE HEATING 


Review of Theory 


The empirical leeward centerline heating method described in reference 4 
Is embodied In the equation 

q . 0.75 q Re 0 ':* (0.002975 + 0.003423 cosA) CD 

c c J£,D 


The parameter q c is the stagnation heating rate on a sphere with a 
radius equal to tfiat of the Shuttle’s upper fuselage. The Reynolds number is 
based on twice that radius. Both quantities are defined by local leeside flow 
properties computed using the flow model shown in figure 2 . It was determined 
that the bow shock angle, 3 » through which free stream flow Is processed 
must be equal to 2a- The flow deflection angle, <5 , depends on and ^ 
3 . The Prandtl-Meyer angle, v , required to expand the flow to the Shuttle s 
upper fuselage is the sum of 6 and ex- 

A pivltol feature of the theory Is the close relationship between axial 
variations in heating and changes In upper fuselage surface flow directions, 

A, measured from oil flew photographs. The technique used to measure oil-flow 
patterns Is illustrated ?n figure 3* The angle, e, between the top 
centerline and a line arawn tangent to the oil-flow path inflection point Is 


* Heating rates reduced from STS— 3 thermocouple data and 3 oj.oi —heating 
corrections for STS— 2 and STS— 3 were provided by D« A. Throckmorton, 

Aero thermodynamics Branch, Space Systems Division, Langley Research Center. 


\ 
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equated with the local angle of attack of flow approaching the top fuselage. 
The local sweep angle is thus defined as 


90° - E 


( 2 ) 


Figure 4 shows two sets of wind-tunnel data where the measured sweep angles 
and the corresponding values of normalized heating rate are plotted as a func- 
tion of x/L at two different test conditions. It is readily apparent that 
sweep angle and top centerline heating vary in an Inverse fashion. A more 
appropriate term for A is "effective'* sweep angle because, as defined for 
the purpose of the top centerline method, its value often becomes larger than 
90°. Effective sweep angles greater than 90° are generally associated with 
inboard flow in the vicinity of the leeward centerline or with surface pat- 
terns caused by flow circulation ahead of the OMS pods and vertical tail. 

This is s impiy a mathematical convenience which allows the theory to penetrate 
zones on the leeward meridian where heating is influenced by various classes 
of sep^raped flow patterns. 4 Using this approach, it was demonstrated that 
the leeward centerline theory is able to cope with the diverse heating envi- 
ronments represented by the wind-tunnel data. As an example, figure 5 shows a 
comparison of the theory's heating predictions with wind-tunnel data at the 
two test conditions indicated in figure 4. 


Extrapolation to Flight Conditions 

The wind-tunnel data base of upper fuselage surface flow directi ons is 
presented In figures 6, 7, and 8 where axial distributions of measured sweep 
angles for Mach numbers of 6, 10, and 14, respectively, are plotted for each 
angle of attack and Reynolds number combination* The hypersonic viscous- 
interaction parameter is also given for each test condition since this term is 
often used to classify the general behavior of leeside separated flow* Sweep 
angles from x/L ■ 0.383 to 0.731 correspond to axial locations where 
thermocouples were positioned on the wind-tunfiel model. Additional measure- 
ments from x/L * 0*30 to 0.787 were made to encompass the locations of orbiuer 
flight instrumentation. For a given Mach number, the Reynolds numbers In 
flight are considerably higher than those in the wind tunnel. Conversely, 
flight Mach numbers are greater than In the wind-tunnel data for corresponding 
Reynolds numbers. The ground-to-f light difference in each of these parameters 
will affect both upper-surface flow patterns and leeward centerline heating. 
Details of the many complex flow Interactions which determine upper fuselage 
heating and surface flow patterns cannot be directly addressed due to the lack 
of information concerning the specific nature of leeside flow processes. 
However, it will be demonstrated that the relatively simple approach described 
here for extrapolating the leeward centerline heating equation and wind-tunnel 
sweep angles to flight conditions is able to capture the essential trends of 
the Reynolds number and Mach number influences on top centerline heating. 

The first step toward converting the empirical leeward centerline heating 
method into a flight prediction technique Is to establish a procedure for 
extrapolating the wind-tunnel sweep angles to their equivalent flight Reynolds 
number values. The next step is to define a criterion which relates the 
flight environment at each trajectory point to the proper set of wind-tunnel 
test conditions In order to duplicate flight trends in leeward centerline 
heating distributions. A third requirement Is to develop a method of 
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►-he wind-tunnel ‘"Mach ^number °f ro^whid^the ^ li tV™™ H “ Ch nUmbers and 

Investigated? 13 *" 6 '* S ° lutl ° n >t0 tbe88 ^ ^ 

may be obtained Sugh C ^l^ulatlor o f h the PPer fUSel3ge 8urface flow patterns 
for each Mach number were cross-^oLrf , swee P~^gle data base. The data 

reveal the format which optimized^! effect of 'Sold “’V* t0 

trends. No easily discernahlp r,r , t6Ct f R£ynolds number related 
found in the Mach 6 data con8istent ^ynolds number trends could be 

trends in the data for a^e fi Till*?™ 1° ^ ° f *-■ ■*»«* 

sweep angle distributions are relativeJ^ind 3 nia J ority of the ^ch 6 

was concluded that the sweep angles meZ - j ependen ^ of Reynolds number. It 
indicated in figure 6 require nf n ured over the Reynolds number range 

Reynolds number? Lweve r f f C °? Ctl0 “ f ° r the effect of 
still be linked to tht IZrZ T f g ° r the different Mach 6 tests must! 
is yet to be described. Similarly ^weenX* 1 conditlon8 in a manner that 

ditions near Jfach 14, where " is ^ / f ° r the tW0 free 8Crean > con- 

for most test cases Til rt ] arge ’ WCra f ° UQd <° be nearly identical 
mediate values of 7* wer! Jubllct lo " directions at «ach 10, for inter- 

of Reynolds number?" Rou!! estimates f variations as a function 

offered on the basis of rh P qo c s of flight Reynolds numbers can be 

>* 1. - o.2 z b o!”‘ tsr ' 7 * ° f , tus 

Mach 6 tests but well below the l n ™ ot m u ,» h maximum value for the 
discussed below, there is a strong Rpvn^f 10 parameter for which, as will be 
the range of 7 *h"I ZZ 0 llT^JT° ldS effect * U PP^ bound on 

perhaps 1.5 oMess. Ihif /f corrections are required, is 
angle distribution for 0 =. 45 ° “V°! tne observation that the axial sweep 
is similar to those at MaH, u t 10 * 16 3nd V " ll219 ln figure 7 

Mach 10 flow pa^eL aS still d°L^t 3ngl r in figure 3. These 

moderate increase in T may dissolve r n n 0n *?“* bUt iC 13 su SS es ted that a 

*• y dlssolv e tne association with Reynolds number. 

plotting sweep Angles atX^Itevnold ^ Reynolds number relationships by 
axial location. ££ £ ^ Ver8U8 * a^ack for each 

fits are used to indicate “ °’J l ia fl ^ re 9 - Unear curve 

each of the three Reynolds numbers t a+°fj th f flow-an 8 le measurements at 
in +10“ of the corresponding^"!’ fit^' £ U8ually fal1 «ith- 

increases the rate of change of the ’ ^ increase in Reynolds number 

attack. The linear curve fits for el^R Mg |j 38 3 functlon of “gle of 

of sweep angle and angle of attack denoted^ 8 numbar share one common value 

flow directions are ELLl 'f by 3nd “r* "*“« 

figure 9 for other axial Iomm Reynolds number. Plots like the one in 

depend on x /l. *^2?“ reV “ - a ^ a fion8 in A r and a , that 

broad sense, A r * r ^ ^ 1C lo! a 

the upper pelage at each axial •utl».^ 1 ^ 1 ££ e ^T P f? lM °" 

the leeside ^ !r° P7 ^ *** ^ ltlaCion of 

££ Si e S2iSS 1 cy lth 

where it begins to ^? rl ^ c ° cjtant unt ii reaching x/L - 0.63 

Patterns with reversed flow directions! pa^uLXy^ 
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and high Reynolds numbers, caused by flow interactions with the blunt forward 
face of the OMS pods and with the vertical tail. The straight-line segments 
in figure 10 which approximate A r are defined as follows: 


for 

0.2 < 

; x/L < 0.34 


(3) 


Ar * 

117.7° - 133.7° 

(x/L) 

for 

0.34 

< x/L < 0.63 




A r * 

72.0° 


(4) 

for 

0.63 

< x/L < 0.8 


(5) 


Ar " 

1*3° + 112.3° 

(x/L) 


shows that ct^- varies linearly and decreases with increasing 
x/L. This is primarily caused by the flow interactions at aft axial stations 
mentioned above. Effective sweep angles within the resulting surface flow 
patterns are elevated at high values of a an< ^ ^co* This forces the poin^ 
of intersection for the various Reynolds number distributions, like those m 
figure 9, to lower angles of attack. The linear curve fit for a r is § ivea 
by 

a r - 40.1° - 29.5° (x/L) (6) 

The parameters A r and a r will be used as part of the procedure to 
correct Mach 10 wind-tunnel sweep angles for the effects of flight Reynolds 
numbers. • 

A basic property of the Reynolds number correction for local flow dirac- 
tions Is related to the slopes of the linear distributions of A vs a- A 
well-defined Reynolds number relationship develops when the slopes of the 
linear curve fits for each Reynolds number, n = AA/Aa, are plotted against 
x/L as in figure 12. The slope of A vs a, n> becomes larger at all axial 
stations with increasing Reynolds number. Each Reynolds number produces a 
maximum value of n at x/L = 0.445 as also shown by the straight-line curve 
fits for the data. The rate of change of n with respect to x/L depends on 
axial location and Reynolds number. This is shown in figure 13 where n' = 
Ati/a(x/L) is presented as a function of log Re*, for locations both fore and 
aft of x/L = 0.445. These terms are defined as follows: 

for x/L < 0.445 

n ' = - 11.059 + 2.542 log Re*, (7) 


for x/L > 0.445 

n ' = 15.429 - 3.297 log Re*, (-) 

The variation of n at x/L = 0.445 provides, the information which allows 
extrapolation of the entire Mach 10 sweep angle data base to higher Reynolds 
numbers corresponding to flight conditions. Figure 14 indicates that the ^ 
three wind-tunnel data points that are available form a linear function of log 
Re*, which is given as 

-13.332 + 2.424 log Re*, 


Hmax = 


(9) 



The dispersion of each data point from the linear correlation is quite scull . 
This would appear to confine the error band of the extrapolation to a very 
narrow range up to Reynolds numbers at least an order of magnitude larger than 
the wind-tunnel values. 

The assessment of Reynolds number effects on sweep angles for the various 
sets of wind-tunnel data and the extrapolation of the Mach 10 upper fuselage 
flow patterns to their equivalent flight Reynolds number form described above 
ccmplete the first step required of applying the empirical leeward-centerline 
heating method to flight conditions. The second step, which was outlined 
earlier , involves the definition of a criterion that relates flight condition, 

to the proper set of wind-tunnel test data. The purpose of this is to "'’sure 

that general trends in heating predictions will reproduce the axial distribu- 
tions of flight heating measurements. The most direct way of obtaining rri s 
information is to collect heating distributions at various entry trajectory 
points and observe which of the wind-tunnel tests produce a corresponding 
inverse variation^ In sweep angles. Examples of such comparisons are illustra- 
ted in figures 15 through 17 for Mach numbers from 24 to approximately 7- 

3oth flights produced nearly the same leeward centerline heating distribu- 

tions at corresponding trajectory points. However, heating rates measurer by 
calorimeters on STS-2 are higher than the thermocouple-derived heat-transfer 
rates for STS-3 at similar flight conditions. This point will be discussed 
later. In each case It Is noted that the correct distribution of sweep tngles 
corresponds to a wind-tunnel test for which the free stream Reynolds number is 
roughly 40 percent of the flight Reynolds number. A study of orbiter leeside 
heating in reference 15 also found that trends in ground-based heating rates 
appeared during entry only at flight Reynolds numbers that were considerably 
higher than An the wind tunnel. Apparently, the wind-tunnel environment more 
accurately simulates flight flow structures on the upper fuselage at higher 
Reynolds numberb over a wide Hach number range. For practical purposes. It is 
sufficient to use the set cf wind-tunnel sweep angles for a freestream 
Reynolds number that is closes;: to 40 percent of the flight value. This 
allows nearly complete coverage of the orbiter's entry trajeptcry instead of 
having heating predictions at only a few discrete flight conditions. This 
relationship seems to be Independent of either wind-tunnel or flight Mach 
number, which indicates that the distribution of leeward centerline heating is 
almost exclusively a function of free stream Reynolds number. According to 
this criterion, top centerline heating predictions at early entry times for 
STS-2 and STS-3 should use wlni-tunnel sweep angles from the Mach 11 data. 
Flight heating distributions from ^ » 20 down to around 10 require the 
extrapolated Mach 10 A f s, Trajectory points below this will use the Mach 6 
wind-tunnel data. 

Heating rates on STS-2 using calorimeters were significantly above the 
STS-3 thermocouple heating measurements. A general dissatisfaction with the 
calorimeters' performance on STS-1 and STS-2 resulted in their removal, ^here 
is also the unsettled question concerning hot-to-cold wall effects on hearing 
rates measured by calorimeters. It can be expected that the large temperature 
differential which existed between the relatively cool calorimeters and the 
surrounding hot surface areas 'would cause these instruments to register a 
higher heating rate than was actually present. The thermocouple data do tor 
suffer from this problem. For these reasons, the determination of a Mach 
number effect on leeward centerline flight heating predictions was based m 
thermocouple measurements from flight 3. 
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Heating predictions based on the Reynolds number extrapolation of wind- 
tunnel sweep angles will be required to obtain the Mach number correction. 

The following is a step-by-step set of instructions on how to ^se true Reynolds 
number correction. Heating predictions for a given trajectory point: will make 
use of sweep angles at the wind-tunnel Reynolds number closest -o 40 percent 
of the flight Reynolds number 

{i0) 

Calculate ntnax and n' for both Reynolds numbers. Next determite n tr - e 
desired axial location using the expressions 

hf = Tw.f - (0.445 - x/Un’f Ola) 

nt " hmax.t " ( 0 . 445 ^- x/L) n ' t (lIb) 

Now compute A for each Reynolds number at the desired angle of attack assuming 
a linear curve fit of A vs a , as in figure 9, using the relations trips 


Af * A r + (a - a r )nf 


(12a) 


A t 53 A r + “ a r^ n t 


(L3) 


The change in sweep angle required to extrapolate the ground-based .il to the 
flight Reynolds number is the difference between Af and A c , thus 


A ext "A + (Af - A t ) Kl-3) 

Where A is interpolated for x/L and a from the win<J-tunnel data sec idett xf ie*- 
by equation (10). Usually, A /= A t because the linear distribution contain- 
ing is only meant to be a general representation of the measurements. 
However, the difference between two such linear representations for a given 
x/L and a is a direct measure of the effect of changing the Reynolds nuaber. 
These steps must be repeated for all axial locations where hearing ^predictions 
are to be made. 


The lack of wind-tunnel data at very high Mach numbers precluded the 
possibility of extracting a Mach number effect on upper fuselage flow direc- 
tions, and the associated leeward centerline heating, from the available 
ground tests. As with the criterion relating flight-heating distributions to 
wind-tunnel surface flow patterns, the effect of flight Mach manber cm the 
extrapolated heating prediction must be formulated using a small puirtion of 
the entry data. Heating rates at = 14.0, * 3.52 x 10 , and =* 

40.9° for flight three were chosen at random and plotted in figure 28(a) along 
with three sets of heating predictions calculated using the MISIVER/*' 0 aero- 
dynamic heating computer program. The highest heating predictions resulted 
from applying the uncorrected wind-tunnel sweep angles for *10*34 and 
Re* =* 1.19 x 10 6 directly to the flight environment. The set pf predictions 
at* intermediate heating levels shows the effect of using the Reynolds number 
extrapolation outlined in equations (10) through (13). The result this 
procedure is a predicted heating distribution that displays the gestural trends 
of the flight measurements, but the predictions are higher by aimosc a factor 
of two. This residual is assumed to be related to the difference between the 
flight and wind-tunnel Mach numbers. It can be accounted for at the Mach 14 
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«;££££*£ .%Krr* heating *— «>• «■»».«. 

y a --ach number correction factor, 4 ,, given as 
i> = (-i» r/’L, r -) 2 

,c ,r (14) 

Ss'red P redI “‘»" s •« good agreement 

consistent with Slight data at ttJ-'t “,T P ™ Ve “ to -' leU results that *rt 
16(b) for JL . 20 0 L - I S , I - .trajectory points, as shown la figure 
wlnd-rannehwe^tgX; ^ “ ; - J J*£-. ponding 

■Jch - ,re plotted as before. There is a auch larged affec- .. 

Kegnojds number ^t”^ ^Kll^ L 

Stt:ri‘ bUtl ° n “ d th ' ^ correction, elation < L) cL Je ‘ 

q c “ °- 73 c >c s 1> (0.002975 + 0.003428 cos A ext ) ( 15 ) 

^6 g or ^ 1 5J^st" ,U J rl ^: ,, L. USe f — P «*!.. ™ fee 

14 tests, A ex (- is assumed to be equal to t \. 


Comparison of leeward Centerline Heating Predictions With Flight Data 
prediction^irh^r 2 ° present comparisons of leeward centerline heating 

he compared . " S^.?2. £U *“ 

d *“ fll8ht 3 ^Vp^.nTto 100 

fppro“nI“"„ ? h " 'll tha l 10 - ltoa “"8 predictions for th. two flight, at 
the slTfe!,, T^i" indl ala"!; “ ndttl,> “ s ™« » S le of attack arc It about 

may be due to Inst rumentafef fectsV rbl/cha”!^ 7 lD geaclEg measurements 
environment Thm - _ , e “ ects r °“- ,,er than large variatiocs In the flirnt 

with 0 ; h magnitude of predicted heating tends to agree more closed 

heat ing ^rate^i^much P the measureme ”^f * ^ P-dicted axial Lstrlbutioi of^ 

for both entry craXJories” 6 £ g measurements of heating distributions 
flight data below Ma C h in on STS-?^^ dl ? a 8 reement between theory ard 

effects, as might be thi case for /*> 20 *“?£ ° f £a strumental 

numbers for StI- 3 are m!ch closer.^ 6 C0, “ PariS0nS a ‘ **«* 

these E l. ld r c " cra " slWo " hrom laminar to turbulent flow can’ be seen In 

^ss*s:;.5'lisirs^ «l \*VZIT‘ da S of n:sht ’• »“ af 
of ** *■ >» £sv~isr r * d 

or laminar flow on the upper fuselage at very high Mach numbers. 
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The theory and flight data rapidly converge for both flight * 

Z ■ 20 . V. entire l-vard -ter « lt ne "Lung lor 

very short time. However, another s y rs no earlier than M* 

STS- 3 concludes that transition to turbulent flow o.curs no ^ wind _ 

- 18. This difference may be related to rhg vlciniry Q f ^ ’ 

tunnel sweep angles used for heating pre - attack low Reynolds "n-xnber 

20. It was noted earlier that tne high ang^ i A 'sweep angles for which 

Mach 10 flow patterns bear a r ® semb ^^ d t0 perhaps by virtue of their 
no Reynolds number correction * 1 te ^ olds n ie: Mach 10 flow angles 

relatively large vaiue of ^ , the i y , -nrrec~ion for 

used in the = 20 heating prediction may T ™ J-tunnei extrapolation. If 
Reynolds number effects than was impose / - <rdica~ a d for the Macn 

so predicted heating rates could be higher than chose cA lea -d * ^ 

20 STS-3 trajectory point. < This would shi w the realm o. t ui 

to somewhat lower Mach numbers But sucn an ^ 

sssrsE.’-srs^- 

antrv heating measurements. 


SIDE FUSELAGE HEATING 


Basic Concepts of Theory 

an empirical method lot predicting aide l(\lt 
the Shuttle orbiter is under development. Li thermocouple heating 

flow patterns and corresponding phase-change \ urbulent 

measurements. The side fuselage t eory u ^ correlation. The 

heating agnation as lor the wind-tunnel f 0 5 % ..«>•. since 

these 1 parameter^are^ close to ' titf ^ciols lot which flight comparisons will 
be made, is 


q 


c 


0.42 



Re 0 ' 3 (0.003531 + 0.004069 cos 
& >D 


A) 


(16) 


from the 


The factor 0.42 correct, the '^o^a’.^To”* lnl« 

stagnation value on a sphere to hpieht of the side fuselage flat 

SSL-M£ SSLM - s- -j - 

“^fSrrre^te’rSarr'taSuo. .10108 only a few data points at 
different values of Re® • 


The surface oil-flow directions 
are also taken to represent angle of 
lage and, thus, a sweep angle in the 


radiating away from the impingement 
attack of flow approaching the siae 
same sense as illustrated in figure 


line 
fuse- 
3 for 



the upper body* The axial variation of sweep angle along the impingement 
location is shown in figure 21 for ^ = 10 and a = 40°. Sweep angles on the 
side fuselage show little change with Reynolds number. They are constant over 
the forward portion of the impingement line and increase sharply at xarge 
values of x/L. The increase in A was determined to be a result of an addi- 
tional expansion of the flow before reaching the fuselage. Both factors 
contribute to a rapid fail in impingement heating at i_hose locations. 

The source of the impinging flow is assumed to be the shear layer which 
originates along a separation line on the strake’s upper surface. It is 
further assumed that separation-point shear-laver flow properties are 
proportional to those at the same axial location ou the st rake ' s leading 
edge. Variations in leading-edge flow properties along the strike were 
accounted for bv interpolation of pressure distributions computed by the 
three-dimensional High Alpha Inviscid Solution (HALIS) 1 ' computer code.* It 
was found that pressure increases linearly along the extent of the strake. 
Another simplifying assumption states that leading edge flow from a given 
fractional distance along the^strake will influence heating at the same 
fractional distance along the side fuselage impingement line. This model 
allows the flow to travel downstream as it moves upward and over toward the 
fuselage . 

These procedures were incorporated into the MINIVER computer program. 
Figure 22 shows an example ox a comparison between the theory’s side fuselage 
impingement line heating predictions and phase-change paint measurements for 
the test condition of =* 10*36, Re„ = 2.37 x 10 b , and a = 40°. The ini- 
tial rise in heating is due to the increase in pressure along the strake £ • * 
leading edge combined with the constant A's in figure 21. Peak heating occurs 
just forward of x/L * 0.4 corresponding to the location at which sweep angles 
begin to increase^. Larger swep angles and the additional expansion of ^ fiow^ 
beyond this point~cause a rapid reduction of impingement heating, Ibe heating 
predictions are in close agreement with wind-tunnel data over the entire 
length of the impingement line. Similar comparisons have been obtained for 
all test conditions at Mach 10 which includes angles of attack from 20° to 40° 
arid free stream Reynolds numbers of 0.59 x 10, 1.19 x 10 , and 2.37 x 10 • 


Comparisons With Flight Data 

The effect of Hach number on the heating prediction has not yet been 
assessed. Therefore, a preliminary comparison with flight data has been 
limited to the STS- 3 trajectory point where * 10.37, Re„ = 5.41 x 10 , 
and a 3 38.9°. Flight Mach number and angle of attack are within the range of 
the wind-tunnel conditions. The flight Reynolds number is larger by over a 
factor of two. But wind-turmel sweep angles, as well as the impingement loca- 
tion, showed little change with Reynolds number at a giver, angle of attack. 
Values of A and impingement location for a s 38.9° were interpolated from the 
wind-tunnel measurements and applied to the flight prediction, figure 23 
shows the resulting axial distribution of impingement heating rate for the 


* Existing HALIS flcv-field computations were supplied by K. James 
Weilmuens ter , Aero thermodynamics Branch, Space Systems Division, Langley 
Research Center. 
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selected flight parameters. The peak, heating rate is approximately ten times 
higher than the average top centerline heating at this same trajectory point. 

It is difficult for the relatively few Instruments on the orbiter's side 
fuselage to obtain a direct measurement of impingement heating. Figure 24 
shows the predicted location of flow Impingement in relation to the positions 
of the side fuselage thermocouples. However, a comparison of the heating pre- 
diction shown in figure 23 with flight data can be accomplished as illustrated 
in figure 25- Cross-sectional measurements cf heating at six axial stations 
are presetted along with the predicted Impingement heating for each value or 
x/L. Wind-tunnel measurements indicate that flow impingement moves off of tr.tr 
upper side fuselage near x/L * 0-65 for a » 40*, so only data for x/L <_ 0 ------ 

are used here. The next downstream array of orbiter instruments Is at x/L = 
0.696. The heating distributions superimposed on the data points were taken 
from wind-tunnel phase-change paint measurements for ^ =* 10.36, Re^ - 
2.37 x l r r at a * 40° and normalized by the predicted Impingement heating. 

The flight heating data generally Conform to the trend and magnitude of the 
projected iistribut Ions . This seems to indicate that the impingement heating 
prediction Is near the correct level. 

Flight lata at x/L « 0.497 and x/L = 0.542 contain some heating measure- 
ments that are approximately 80 and 120 percent above the mean loca] values , 
respectively. Figure 26 shows that these pulses of high convective heating 
are associated with large and erratic excursions In surface temperature which 
occur at later entry times. Similar temperature fluctuations affect many of 
the side fuselage thermocouple locations at slightly different times, but 
these variations are always confined to the high Reynolds number portion of 
the trajectory. The anomaly at x/L * 0.497 is above the impingement line 
while at x/L « 0.594 it is well below the impingeiaent location. Many similar 
heating "spikes" were observed in the wind-tunnel data for Mach 6 and Mach 10, 
but only at locations above the impingement line. This is Illustrated in 
figure 27 using cne of the Mach 6 test cases for which Reynolds number is the 
same as for the STS-3 trajectory point. The average increase over local 
heating associated with these features in the wind tunnel was about 70 
percent. The phase-change paint data revealed that the phenomena are highly 
localized, ar indicated by the slender heating profile In figure 27. This 
same profile was dpplied to the data In figure 25 using dashed lines In order 
to distinguish actual flight measurements from normalized wind-tunnel heating. 

It Is suggested that these elevated local heating rates are caused by 
embedded vortices which are generated by viscous interactions during the 
impingement process. Embedded vortices are believed to be caused by boundary- 
layer cross-flow Instabilities. 18 References 19 and 20 are examples of 

the many studies on the relation between embedded streamvise vorticity and 
flow Impingement. Figure 28 contains a photograph showing a sequence of 
uniformly spaced streaks in phase change paint above, and originating from, 
the imping ement location on the orbiter model's side fuselage at a Mach 6 test 
condition. Each streak is thought to represent a very thin line of vortex 
impingement which produces locally higher heating, and the !uagnltude of heat- 
ing decreases along its length. Side fuselage "streak" heating has also been 
noced on an early phase B orbiter conf iguratioj^ ' and on the ASSET entry 
vehicle. 22 A larger number of streaks were observed in Che pnase-change paint 
tests with increasing Reynolds number and higher angles cf attack. Streaks 
vere present for all test conditions except a * 20° and 30* for Mach 10. 
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_ 6Se re ” 3 aai cate that side fuselage "streak" heating nay well be expected 
occur in the entry flight regime containing the selected STS-3 trajectory 
point. In addition, the changing number of streaks and variations in the 
spac ng between them at different free-stream conditions and angles of attack 
mean that embedded vortices on the orbiter in flight will move longitudinally 
on the side fuselage. This motion will cause a number of individual vortices 
to sweep across a fixed location resulting in intermittent locally higher 
. a ‘ could explain the large temperature variations at later entry 

imes shown in figure 26. Very large side fuselage STS-2 heating rates that 

vorticity prevlously docu nented may also be caused by the onset of embedded 


CONCLUDING REMARKS 


A method has been developed for extrapolating a wind-tunnel-developed 
empirical heating technique for the Space Shuttle orbiter' s leeward centerline 
° light conditions. The distribution of heating along the vehicle's leeward 
meridian was found to be primarily a function of Reynolds number. Axial heat- 
ing trends in t light correspond to those In the wind tunnel for which the 
Reynolds number is approximately 40 percent of the flight value. Only those 
ind tunnel leeside tuselage flow patterns at intermediate values of 7 
displayed significant and consistent sensitivity to changes in Reynold! 
number. The effect of Mach number on heating predictions was determined 

roug mited use of the flight data. Application of the extrapolated heat- 
ing method to the flight environments of STS -2 and STS-3 produced generally 
favorable comparisons. It was tentatively concluded that the STS-2 heatin- 
measurements were of lower quality than those of STS-3. The theory may" *** 
provide a somewhat conservative indication for the time of transition from 

t0 e tu . r j> ulent flow - Heating predictions afforded by this procedure are 
adequate for the design of upper fuselage thermal protection systems. 

* n ® w technique for computing side fuselage Impingement heating was 

briefly described. This method is derived from the leeward centerline 

t eory. Although still under development, the side fuselage heating method 
was s own to z^ree well with wind-tunnel data and with selected STS-3 flight 
measurements. The comparison with flight data revealed that, as in the wind 
tunnel, there are areas of locally enhanced heating at side fuselage locations 
well away from the Impingement line, flie associated heating rates were 
approximately 100 percent higher than nearby undisturbed levels. It is 
suggested that this phenomenon Is caused by embedded vortices resulting from 
viscous interactions that are perhaps related to flow reattachment at 
tree-stream conditions which satisfy critical values of Mach number and 
Reynolds number at a given angle of attack. The existence of these features 
will have an impact on thermal protection requirements of future winged entry 
vehicles which experience flow impingement cn the side fuselage* * 
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field model used to compute local leeside properties 
or leeward centerline heating theory. 
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Figure 3.- Upper fuselage surface pattern showing measurement of flow angle. 



Figure 4.- Illustration of the inverse relation between effective sweep angle 
and leeward centerline heating rate. 
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Figure 9.- Illustration of Reynolds number effect 
on Mach 10 wind tunnel sweep angles. 
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Figure 15.- Comparison of flight leeward centerline heating dis t rLbu t ions 
near = 24 with wind tunnel ef r ective sweep angle variations. 
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Figure 16. Comparison of flight leeward centerline heating distributions 
near = 14 with wind tunnel effective sweep angle variations- 
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Figure 17.- Comparison of flight leeward centerline heating distribu 
near M = 7 with wind tunnel effective sweep angle variations. 
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Figure 18.- Comparison of flight leeward centerline heating rates 
with predictions in various stages of correction for the 
effects of Mach number and Reynolds number. 
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Figure l^.- Comparison of leeward centerline heating predictions 
with STS-2 flight data. 
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Figure 21.- Effective sweep angles resulting 
from side fuselage impingement. 
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Figure 22.- Example of comparison between side fuselage impingement 
heating predictions and wind tunnel measurements for = 10.36, 
Re = 2.37 x IQ 6 and a = 40°. 
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Figure 23.- Predicted side fuselage impingement heating for STS-3 
at M = 10.37, Re = 5.4 x 10 6 , and a = 38.9°. 
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Figure 24.- Predicted location of impingement line with respect to side 
fuselage thermocouples at M = 10.37 and a - 38.9° during STS-3. 
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(a) x/L = 0.497; s/L = 0.0950; V07T9924. 



(b) x/L = 0.542; s/L = 0.1248; V07T9905. 

Figure 26.- STS-3 side fuselage temperature time histories shoving 
temperature fluctuations for entry times near M = 10.37. 
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Figure 27.- Example of heating spike 
M = 6, Re = 5.4 * 10 6 , and a = 
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in wind tunnel data for 
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